Hawking-like radiation does not require a trapped region 



O 
O 

(N 



C. Barcelo, 1 S. Liberati, 2 S. Sonego, 3 and M. Visser 4 

1 Institute) de Astrofisica de Andalucia, CSIC, Camino Bajo de Huetor 50, 18008 Granada, Spain 
2 International School for Advanced Studies, Via Beirut 2-4, 
34014 Trieste, Italy and INFN, Sezione di Trieste 
3 Universita di Udine, Via delle Scienze 208, 33100 Udine, Italy 
^School of Mathematics, Statistics, and Computer Science, Victoria University of Wellington, New Zealand 

(Dated: 22 June 2006; WI^K-ed February 4, 2008) 

We discuss the issue of quasi-particle production by "analogue black holes" with particular at- 
tention to the possibility of reproducing Hawking radiation in a laboratory. By constructing simple 
geometric acoustic models, we obtain a somewhat unexpected result: We show that in order to 
obtain a stationary and Planckian emission of quasi-particles, it is not necessary to create a trapped 
region in the acoustic spacetime (corresponding to a supersonic regime in the fluid flow). It is 
sufficient to set up a dynamically changing flow asymptotically approaching a sonic regime with 
sufficient rapidity in laboratory time. 
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Introduction : It is by now well established that the 
physics associated with classical and quantum fields in 
curved spacetimes can be reproduced, within certain ap- 
proximations, in a variety of different physical systems 
by constructing so-called "analogue spacetimes" Q, Q. 
Among such systems the simplest example is provided by 
a barotropic, irrotational and viscosity-free fluid. Small 
classical (or quantum) disturbances propagating in this 
type of fluid behave equivalently to a linear classical (or 
quantum) field over an effective acoustic spacetime, en- 
dowed with an acoustic metric. In the case of a fluid flow 
along an infinitely long thin pipe, with density and ve- 
locity fields constant on any cross section orthogonal to 
the pipe, one can write the (1 + l)-dimensional acoustic 
metric as 

g= Q 2 [- (c 2 -v 2 )dt 2 -2vdtdx + dx 2 ] , (1) 

where t € H denotes laboratory time and x € H the 
physical distance along the pipe. In this expression c is 
the speed of sound, v is the fluid velocity, and SI is an un- 
specified non- vanishing function U . In general, all these 
quantities depend on the laboratory coordinates x and t. 
Here, we shall assume that c is a constant. Hence, it is 
the velocity v(x, t) that contains all the relevant informa- 
tion about the causal structure of the acoustic spacetime. 
(See Ref. Q| for a detailed analysis of the causal struc- 
ture associated with a broad class of (1 + l)-dimcnsional 
acoustic geometries, both static and dynamic.) 

In this letter we analyze, in simple terms, the is- 
sue of quantum quasi-particle creation (in this specific 
case, phonon creation) by externally specified (1 + 1)- 
dimensional analogue geometries simulating the forma- 
tion of black hole-like configurations. (Several additional 
cases and a more extensive discussion of the results can 
be found in Ref. Q.) In this analysis we have in mind, 
on the one hand, the possibility of setting up laboratory 



experiments exhibiting Hawking-like radiation 0, Q and, 
on the other hand, the acquisition of new insights into the 
physics of black hole evaporation in semiclassical gravity. 
In particular, we have found, quite surprisingly, that in 
order to produce a Hawking-like effect it is neither nec- 
essary to generate a supersonic regime (fluid velocity v 
strictly larger than sound velocity c), nor even a sonic 
point at finite laboratory time. All one needs is that a 
sonic point develops in the asymptotic future (that is, for 
t — » +00) with sufficient rapidity (we shall in due course 
explain exactly what we mean by this). 
General framework: A sonic point in the flow, where 
v(t, x) = ±c, corresponds to a so-called acoustic trapping 
(or apparent) horizon for the Lorentzian geometry de- 
fined by the metric (JTJ. See e.g. Ref. Q, Sec. 2.5, pp. 15- 
16. Consider a monotonically non-decreasing function 
v(x) such that v(0) = —c and v(x) — > for x — > +00. 
If one chooses v(x, t) = v(x) in Q, the corresponding 
acoustic spacetime represents, for observers with x > 0, 
a static black hole with the horizon located at x = (in 
this case trapping and event horizon coincide), a black 
hole region for x < 0, and a (right-sided) surface grav- 
ity i 



:= lim 



dv(x) 
dx 



(2) 



Now, taking the above v(x) 
dependent velocity functions 



let us consider t- 



u(x t) i if X ^ m ' 



(3) 



with £ a monotonically decreasing function of t, such 
that lim £(t) = +00 and lim £(t) = 0. (The first 

t — > — oc t — > — 00 

condition serves to guarantee that spacetime is flat at 
early times, whereas we impose the second one only for 



FIG. 1: The plot shows a series of temporal snapshots, ti < t2 < tz < ti, of the velocity profile v(x,t). In the first-row series, 
the sonic regime is reached at a finite laboratory time (this series represents the formation of a non-extremal acoustic black 
hole). In the second-row series, the sonic regime is reached only asymptotically in laboratory time (this series represents the 
formation of what we called a critical acoustic black hole). 



simplicity.) The derivative with respect to x of the cor- 
responding velocity profile is discontinuous at x = 
where a "kink" traveling to the left is located (see Fig.^l. 
There are basically two possibilities for £, according to 
whether the value £ = is attained for a finite laboratory 
time £h or asymptotically for an infinite future value of 
laboratory time. 

In the first case £(£h) = and the corresponding metric 
represents the formation of a non-extremal black hole 
(in this letter we will consider only the case k ^ 0, ex- 
tremal black holes — zero surface gravity — are treated 
m @)- F° r small values of \t — %| we have 



£(t) = -X(t-tn) + 0([t-t n } 



(4) 



where A is a positive parameter. Apart from this feature, 
the detailed behaviour of £ is largely irrelevant for our 
purposes. 

If instead £ — » is attained only at infinite future time, 
that is lim £(t) = 0, one is describing the asymptotic 

t — >+oc 

formation of what we shall call a critical black hole. In 
the following we shall restrict the discussion to the case in 
which for t — > +oo, the behaviour of £ is exponential, that 
is £(t) ~ Ae~ Knt , with kd a positive constant, in general 
different from k, and A > 0. Of course, other possibilities 
for the asymptotics can be envisaged; for example, in 
Ref. []| the case of a power law £(t) ~ B t~ v , with v > 
and B > 0, is also studied in detail. 

Starting with a quantum scalar field in its natural 
Minkowskian vacuum at t — > — oo, we want to know the 
total quantity of quasi-particle production to be detected 
at the right asymptotic region at late times, t — ► +oo, 
caused by the dynamical evolution of the velocity pro- 
file v(x,t). This can be done by following the standard 
procedure in which one computes the Bogoliubov (3 coef- 
ficients between modes that have positive frequency with 
respect to null coordinates U and u, with U regular on 



the acoustic event horizon TL, whereas u tends to +oo as 
H is approached. (See, e.g., Ref. We must thus find 
the relation between U and u for a sound ray that is close 
to the horizon, i.e., in the asymptotic regime u — > +oo. 
If this relation turns out to be exponential, it is then a 
well established result that a Planckian spectrum will be 
observed at late times 0, so Hawking-like radiation is 
recovered. 

In the geometric acoustic approximation, a right-going 
sound ray is an integral curve of the differential equation 



d:r 



c + v(x, t) , 



(5) 



which can be integrated exactly in order to find the values 
of U and u for a generic right-moving ray that crosses, 
at laboratory time to, the kink x — £(t) (details can be 
found in Ref. 0): 



U = tn- 



d* «(£(*)) 



tn - 



£(*d) 



c Jew c 



v(x) 



(6) 



(7) 



Eliminating to between Eqs. © and provides the 

sought-for relation between U and u. 

Formation of a non- extremal black hole: In the case of 

the formation of a non-extremal black hole, the function 

v(x) must be such that for small values of \x\, one can 

write 



v(x) = —c + Kx + 0(x 2 ) 



(8) 



In addition, the trapping horizon forms at a finite labo- 
ratory time, say at t = tn, just when £(4h) = 0. In this 
situation an event horizon always exists, generated by 
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the right-moving ray that eventually remains frozen on 
the trapping horizon, at x = 0. For such a ray to — > tn, 
and since = 0, the U parameter has the finite value 



Un = tn + \ I (1/ '•((.(/)) 



(9) 



For a ray with U < C/h we then obtain, combining 
Eqs. © and ©: 



U = U B + t - t H - 



tH 



dttJ(£(t)) • (10) 



Co 



This ea;act equation is now in a form suitable for con- 
veniently extracting approximate results in the region 
~ ^h, corresponding to sound rays that "skim" the 
horizon. Using Eqs. (QJ and © we find U = J7 H + 7 (*o - 
tu)+0([to—tn] 2 ). This provides us with the link between 
U and to- 

In order to link to with u, consider the integral on 
the right hand side of Eq. J7J. For x —> +00, the inte- 
grand function vanishes, while near £(to) it can be ap- 
proximated by —c/(kx). Then the integral is just given 
by the difference of the corresponding integrals evaluated 
at x = +00 and x — £(to), respectively, up to a possi- 
ble finite positive constant. This gives — A (to — £h) ~ 
const e~ KlU . Together with the previous link between U 
and to this leads to 



U ~ ?7h — const e 



(11) 



This relation between U and u is exactly the one found by 
Hawking in his famous analysis of particle creation by a 
collapsing star /6| . It is by now a standard result that this 
relation implies the stationary creation of particles with 
a Planckian spectrum at temperature T# = n/(2ir) 
Formation of a critical black hole: Consider now the 
same type of function v(x) but this time the sonic point is 
approached asymptotically in an infinite amount of labo- 
ratory time. Now the trapping horizon is just an asymp- 
totic point located at x — 0, t — > +00, and in order to 
establish whether an event horizon does, or does not, ex- 
ist one must perform an actual calculation of Uu for the 
"last" ray that crosses the kink. The expression for Uu is 
again obtained from Eq. |JB}, noticing that now to = +°° 
along the generator of the would-be horizon, so 



lim ( to 

t ->+oo 



dtvm) 



(12) 



The necessary and sufficient condition for the event hori- 
zon to exist is that the limit on right hand side of Eq. I|12|) 
be finite. The integrand on right hand side of (|12f> can 
be approximated, for t — > to — > +00, as — c + «£(t), while 
for t — > — 00 it just approaches zero. Hence Uu is, up 
to a finite constant, equal to n/c times the integral of £, 
evaluated at t — * +00. For an exponential behaviour, Uu 
turns out to be finite. 



For another right-moving sound ray that corresponds 
to a value U < Uu, combining Eqs. © and flT2")l. realizing 
that v ~ — c + n£(t) for £(t) close to zero, and using 
expansion |JSj, we find 

[/^ H _iM_* / + °°dt£(t). (13) 



For the asymptotically exponential £ (t) considered in this 
letter this implies 



For the link between to and u we obtain 
u ~ t - - ln£(i ) , 



(14) 



(15) 



as one can easily check by inserting the appropriate 
asymptotic expansions into Eq. J7J). Using Eqs. i|15|) and 
l[rl)l we finally find 



U ~ Uu — const exp u 

k + K D 



(16) 



Therefore, we have shown that by setting up a fluid flow 
that reaches a sonic regime asymptotically with expo- 
nential rapidity in laboratory time one can produce a 
stationary and Planckian spectrum of quasi-particles at 
asymptotic infinity with temperature T c ff = K e ff/(27r), 
where 



«eff 



K KD 



(17) 



When kd ^> k, « c ff ~ k and thus, our result is op- 
erationally indistinguishable from that of Hawking. In 
this sense neither the existence of an ergoregion nor that 
of an apparent horizon are needed for the simulation of 
Hawking radiation. Remarkably this is true also for a 
"double-sided" velocity profile in which the asymptotic 
sonic regime is limited to a single spatial point Q . 

In terms of inverse temperatures, defined as /3, — 
1/Tj = 2ir/ '«;,-, we have the very suggestive result 



Pcft = /^Hawking + 0D ■ 



(18) 



Experimental realizability: The critical black hole model 
seems worth taking into consideration in connection with 
the realizability of a Hawking-like flux in the laboratory. 
The creation of supersonic configurations in a laboratory 
is usually associated with the development of instabili- 
ties. There are many examples of the latter in the liter- 
ature; e.g. in Ref. |9j it was shown that in an analogue 
spacetime based on ripplons on the interface between two 
different sliding superfluids (for instance, 3 He-phase A 
and 3 He-phase B), the formation of an ergoregion would 
make the ripplons acquire an amplification factor that 
eventually would destroy the configuration. Therefore, 
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this analogue system, although very interesting in its 
own right, will prove to be useless in terms of detect- 
ing a Hawking-like flux. However, by creating, instead of 
an ergoregion, a critical configuration one should be able 
to at least have a better control of the incipient insta- 
bility, while at the same time producing a dynamically 
controllable Hawking-like flux. 

Nevertheless, the actual realization of a critical con- 
figuration might also appear problematic for entirely dif- 
ferent reasons. The corresponding velocity profiles are 
characterized by discontinuities in the derivatives, so one 
might wonder whether they would be amenable to exper- 
imental construction, given that the continuum model is 
only an approximation. In particular, could it be that in 
a real experimental setting one does not observe the re- 
sults predicted by our analysis, but those corresponding 
to models where the discontinuities have somehow been 
smoothed out? 

Fortunately, this will not be the case. In realistic situ- 
ations, what is relevant for Hawking radiation is a coarse- 
grained profile — obtained by averaging over a scale A 
larger than the one which characterizes the breakdown 
of the continuum model — which does not contain the 
unphysical small scale details of v(x, t) . This implies that 
the reliable results are those involving features recogniz- 
able over length scales of order A, or larger [3- In par- 
ticular, the relevant surface gravity will be defined by av- 
eraging the slope of the velocity profile over scales which 
are of order of A. This averaged surface gravity will be 
non-zero for the critical profile, as well as for a smoothed 
one. Indeed, it will be approximately equal to the surface 
gravity at the horizon of the critical black hole, while it 
will obviously not coincide with the one of the smoothed 
profile (which is zero) [jUl^]. 

Hints for semiclassical gravity: The critical collapse re- 
sult also suggests an alternative scenario for the semiclas- 
sical collapse and evaporation of black hole-like objects. 
Although somewhat speculative at this stage we believe 
it is worthwhile to describe it here. 

Imagine a dynamically collapsing star. The collapse 
process starts to create particles dynamically before the 
surface of the star crosses its Schwarzschild radius. (This 
particle creation is normally associated with a transient 
regime and have nothing to do with Hawking's Planck- 
ian radiation.) The energy extracted from the star in 
this way will (due to energy conservation) reduce its to- 
tal mass, and so also its Schwarzschild radius, so k in- 
creases. By this argument alone, we can see that a pro- 
cess is established in which the surface of the star starts to 
"chase" its Schwarzschild radius while both collapse to- 
wards zero (this situation was already described by Boul- 
ware in Ref. [jof). 

Now, contrary to the standard view, it could happen 
that the back-reaction on the geometry is such that it 
prevents the surface of the star from actually crossing its 
Schwarzschild radius (see 0] for a discussion of this pos- 



sibility). The function £(i) in our calculation would rep- 
resent the radial distance between the surface of the star 
and its Schwarzschild radius. Then it is sensible to think 
that during the evaporation process any dynamical kd 
would also depend on t. As the evaporation temperature 
increases (k increases) the back-reaction would become 
more efficient and therefore we might expect that Kr> de- 
creases. Then the evolution of the evaporation temper- 
ature would interpolate between an early-time tempera- 
ture completely controlled by n(t), and a late-time tem- 
perature completely controlled by «d(^)- Therefore, this 
black hole-like object will have a temperature-decreasing 
phase, showing a possible semiclassical mechanism for 
regularizing the end point of the evaporation process. 

In this scenario the complete semiclassical geometry 
will have neither a trapped region (and so by construc- 
tion no trapping horizon) nor an event horizon. In this 
circumstance there would be no trans-Planckian prob- 
lem, nor information loss associated with the collapse 
and evaporation of this black hole-like object. Whether 
this scenario is viable or not will be the subject of future 
work. 
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